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Abstract 

This paper answers a question of Emery (2009) by constructing an explicit coupling of two copies of the Be- 
nes, Karatzas, and Rishel (1991) diffusion (BKR diffusion), neither of which starts at the origin, and whose 
natural nitrations agree. The paper commences with a brief survey of probabilistic coupling, introducing 
the formal definition of an immersed coupling (the natural filtration of each component is immersed in a 
common underlying filtration; such couplings have been described as co-adapted or Markovian in older 
terminologies) and of an equi-filtration coupling (the natural filtration of each component is immersed in 
the filtration of the other; consequently the underlying filtration is simultaneously the natural filtration 
for each of the two coupled processes). This survey is followed by a detailed case-study of the simpler but 
potentially thematic problem of coupling Brownian motion together with its local time at 0. This prob- 
lem possesses its own intrinsic interest as well as being closely related to the BKR coupling construction. 
Attention focusses on a simple and natural immersed (co-adapted) coupling, namely the reflection / syn- 
chronized coupling. It is shown that this coupling is optimal amongst all immersed couplings of Brownian 
motion together with its local time at 0, in the sense of maximizing the coupling probability at all possible 
times, at least when not started at pairs of initial points lying in a certain singular set. However numerical 
evidence indicates that the coupling is not a maximal coupling, and is a simple but non-trivial instance 
for which this distinction occurs. It is shown how the reflection / synchronized coupling can be converted 
into a successful equi-filtration coupling, by modifying the coupling using a deterministic time-delay and 
then by concatenating an infinite sequence of such modified couplings. The construction of an explicit 
equi-filtration coupling of two copies of the BKR diffusion follows by a direct generalization, although the 
proof of success for the BKR coupling requires somewhat more analysis than in the local time case. 
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1 Introduction 

We begin with a brief survey of probabilistic coupling, which serves both to introduce some key concepts and 
to establish a context for the results proved in this paper. The concept of coupling has a long and distinguished 
history, dating back to Doeblin (1938) (a biographical appreciation is given by Lindvall, 1991). The method is 
now the subject of two scholarly expositions (Lindvall, 1992; Thorisson, 2000), and has become a standard tool 
of the working probabilist. Historically the thematic problem for coupling is that of constructing two copies of 
a given process on the same sample space, starting at two different starting points but eventually coinciding. 
Such a coupling is said to be successful. In fact many applications of coupling avoid the notion of eventually 
coinciding (this more general concept appears in ergodic theory as the notion of a "joining"); nevertheless 
the thematic problem has been formative for the theory and remains significant in developing methods and 
intuition. Probabilistic coupling in general has found application throughout probability, for example in 
construction of gradient estimates, in distributional approximation (for instance Stein-Chen approximation), 



in perfect simulation, and in monotonicity results for heat equations in insulated domains. The study of 
coupling in its own right is therefore a foundational topic for probability theory. 

A landmark development in the study of coupling was the introduction of the notion of maximal coupling: 
a coupling which simultaneously maximises the chances of succeeding before time t for all possible t. Perhaps 
it will surprise the reader to learn that maximal couplings always exist: this was established by Griffeath (1975) 
for time-homogeneous discrete Markov chains and by Goldstein (1978) for more general discrete-time processes, 
based on a tail cr-algebra condition. (Note that even a maximal coupling need not necessarily have probability 
1 of succeeding!) See also the very explicit construction given by Pitman (1976) for time-homogeneous discrete 
Markov chains, Sverchkov and Smirnov (1990)'s note on coupling for continuous time using the Ji topology, and 
Thorisson (1994)'s notion of shift coupling, which weakens the coupling requirement by allowing for general 
time-shifting of the coupled processes. (An informative treatment of some of the subtleties is given in the 
treatment of "faithful coupling" in Rosenthal, 1997.) 

In general the construction of maximal couplings is a demanding business: for substantial applications the 
task of construction is liable to require at least as much knowledge of the process in question as might be 
needed to solve the original problem to which the coupling method is to be applied. (Notwithstanding this 
general and justifiable pessimism, the simple reflection coupling of Brownian motion is a successful maximal 
coupling. Attention was originally drawn to this construction in the very influential unpublished preprint 
of Lindvall, 1982.) More commonly one works with more constructable but less powerful couplings, such as 
"co-adapted couplings". Co-adapted couplings (sometimes also called "Markovian couplings" in the context of 
coupling of Markov processes) require the two copies of the processes concerned to be adapted to the same 
filtration, and to have the same conditional laws based on conditioning on filtration cr-algebras. In the succinct 
language of nitrations (c.f. Beghdadi-Sakrani and Emery, 1999; Emery, 2005, 2009), the natural nitrations 
of the two processes must both be immersed in a common filtration (that is, the martingales of the natural 
nitrations must remain martingales in the larger common filtration). We therefore propose and adopt the new 
terminology of immersed couplings to replace the nomenclature of co-adapted or Markovian couplings: 

Definition 1. Consider two processes X and Y. An immersed coupling of X and Y is a construction of copies X, 
Y of X, Y, defined on the same probability space (CI, 5", P), and adapted to the same filtration {SF t : t ^ 0}, such 
that any martingale in the natural filtration of X remains a martingale in the common filtration {Jt : t ^ 0}, 
and likewise for any martingale in the natural filtration of Y. 

The extent to which immersed couplings are less powerful than maximal couplings finds an initial assessment 
in Burdzy and Kendall (2000), where they are studied in the guise of Markovian couplings. As part of a study 
of shy coupling (the antithesis of the thematic coupling problem, in which one seeks to construct coupled 
copies which almost surely stay at least a fixed positive distance apart), Kendall (2009, Lemma 6) records a 
characterization of immersed couplings of Brownian motion which has long been part of the general folklore of 
stochastic calculus: any immersive coupling of two d-dimensional Brownian motions A and B can be represented 
by the stochastic differential equation 

dA = J T dB + K T dC (1) 

where C is a Brownian motion independent of B (perhaps to be defined after augmenting the filtration, if this 
is necessary to construct C), and J and K are two (dx d) matrix- valued predictable random processes satisfying 
J T J + K T K = I where I is the (d x d) identity matrix. We can view J as a predictable matrix- valued control for 
a somewhat degenerate stochastic control problem. (An informal discussion of links between stochastic control 
and coupling can be found in Kendall, 2007, Section 2.) 

The terminology of immersed couplings is useful not only for its succinct definition, but also because it 
draws attention to a stricter constraint. Additionally one could demand that either of the coupled copies 
could be constructed from the other, which corresponds to the requirement that the coupling possesses the 
equi-filtration property: 

Definition 2. Consider two processes X and Y. An equi-filtration coupling of X and Y is an immersed 
coupling X, Y such that the natural filtration of X is immersed in that of Y, and vice versa. 
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Of course it is the case that the equi-filtration coupling property follows from each natural filtration being 
immersed in the other. Consider the simplest nontrivial example of diffusion coupling; Lindvall's Brownian 
reflection coupling is not only a successful maximal coupling, but also immersed and even equi-filtration. This 
is a very special case; for example Connor (2007, PhD. thesis) notes that reflection coupling of the driving 
Brownian motions is clearly immersed but is not maximal even in the simple case of the Ornstein-Uhlenbeck 
process, if one copy of the Ornstein-Uhlenbeck process is started from and the other copy is started from 
equilibrium. (Further exploration of the difference between maximality and immersion for couplings can be 
found in Kuwada and Sturm, 2007; Kuwada, 2009.) 

The first objective of this paper is to investigate and explore properties of the construction of immersed 
and equi-filtration couplings in the simple case of coupling Brownian motion together with local time at 0. 
As a coupling problem this is only a little more complicated than the basic Brownian motion case, but it 
produces an example of existence of a successful immersed coupling (the reflection / synchronized coupling, 
Definition 3) which is strictly optimal among all immersed couplings but (according to numerical evidence) is 
not maximal (Theorems 9, 10, 13 below). The reader may wish to compare other work on optimal immersed 
couplings for random walks on the line, on hypercubes and on hypercomplete graphs (Rogers, 1999; Connor 
and Jacka, 2008; Connor, 2009). 

A significant motivation for this study arises from the consideration that the reflection coupling has been 
a model for a wide variety of more sophisticated immersed couplings. For example, reflection coupling has 
been generalized to the case of elliptic diffusions with smooth coefficients (Lindvall and Rogers, 1986; Chen 
and Li, 1989), and also to the case of Riemannian Brownian motion (Kendall, 1998), in which case there are 
connections with curvature properties. More recently coupling techniques have been extended to cover some 
cases of hypoelliptic diffusions (Ben Arous, Cranston, and Kendall, 1995; Kendall and Price, 2004; Kendall, 
2007, 2010); essentially the issue here is to couple simultaneously not only Brownian motion but also one 
or more path functionals of the Brownian motion, namely time integrals, iterated time integrals, and Ito 
stochastic area integrals. Here it is necessary to augment the reflection coupling strategy with other coupling 
strategies, notably synchronous coupling and rotation coupling. In the stochastic differential framework (1), 
synchronous coupling corresponds to K — and J = I, while rotation coupling corresponds to K = and J 
equal to a d-dimensional rotation. (It is interesting to compare this direction of research with the work of 
Emery, 2005, Theorem 1; this characterizes Brownian nitrations using the notion of "self- coupling" - jointly 
immersed Brownian nitrations for which a prescribed scalar functional is approximately coupled.) 

While Brownian motion together with local time at does not form a hypoelliptic diffusion in the strict 
sense, nevertheless the question of its coupling theory is clearly related to the hypoelliptic couplings mentioned 
above. The successful reflection / synchronized coupling is not only simple, but also (in view of the results 
proved here) evidently the right coupling for this situation. It is reasonable to hope that a careful and complete 
study of the reflection / synchronized coupling will be helpful in formulating and studying coupling methods 
for more general situations, as well as suggestive for coupling theory for hypoelliptic diffusions. The second 
objective of the paper is to demonstrate the first fruits of this aspiration and is fulfilled in Theorem 19 below, 
exhibiting a successful equi-filtration coupling for the BKR diffusion. The approach follows closely the methods 
developed for the reflection / synchronized coupling for Brownian motion together with local time at 0. 

In summary, then, this paper conducts a case study of an almost surely successful coupling of a simple non- 
elliptic diffusion in the context of immersed and equi-filtration couplings; namely the reflection / synchronized 
coupling for Brownian motion together with local time at 0. The results of this case study are then applied 
to answer a question raised by Emery (2009), by constructing an explicit equi-filtration coupling for BKR 
diffusions neither of which are begun at the origin. 

Section 2 introduces the simple reflection / synchronized coupling for Brownian motion together with local 
time at 0, exploiting Tanaka's formula and the Levy transform to re-cast the problem in terms of coupling 
Brownian motion together with its running supremum. The simplicity of this coupling allows for explicit 
calculation: in particular it is shown that the reflection / synchronized coupling is strictly optimal amongst all 
immersed couplings, at least when their starting conditions are non-singular (here "optimal" means optimal in 
the sense of maximizing the probability of coupling by a given time t, for all possible times t, while "singular" 
means that in the re-cast form the two running suprema processes do not start from the same level). The 
moment-generating function for the coupling time is computed, and compared numerically with the moment- 
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generating function for the maximal coupling time: numerical calculation then indicates that the reflection / 
synchronized coupling cannot be a maximal coupling. 

The reflection / synchronized coupling is an immersed coupling but is not equi-filtration. Section 3 shows 
that if the couplings are perturbed by a simple deterministic time delay then it is possible to use a sequence 
of the resulting approximate couplings to construct a successful equi-filtration coupling of Brownian motion 
together with its local time at 0. 

Section 4 introduces the BKR diffusion, sketches the immersed coupling described in Emery (2009) (which 
bears a strong family resemblance to the reflection / synchronized coupling of Section 2, and which therefore 
is described here as a variant reflection / synchronized coupling), and notes that significant components of 
this variant reflection / synchronized coupling are actually immersed in the natural nitrations of both coupled 
diffusions. This is used to generate a successful equi-filtration coupling using the strategy of Section 3, hence 
answering Emery's question. 

The paper is concluded by Section 5, which reviews the results of the paper and discusses some further 
research questions. 

2 Coupling Brownian motion together with local time 

The purpose of this section is to exhibit a successful immersed (but not equi-filtration) coupling for the two- 
dimensional diffusion made up of Brownian motion together with local time at zero. The simple construction 
(known already to Emery, 2009) is based on Tanaka's formula for Brownian local time, and permits informative 
exact computations. In particular we are able to prove strict optimality of this coupling amongst all immersed 
couplings (Theorem 9), so long as the initial conditions are non-singular in a manner to be explained below, 
and thus to establish the optimal rate of immersed coupling (Theorem 10). 

2.1 Representation via the Tanaka formula 

Recall the Tanaka formula or Levy transform, expressing Brownian local time at in terms of a stochastic 
integral: 

d|X| = sgn(X)dX + dL (0) . (2) 

Here X is a real Brownian motion and L' ^ is the local time accumulated by X at 0. An immediate consequence 
of (2) is Levy's famous transform, which represents |X| and L^ ' in terms of a new real Brownian motion B and 
its running supremum S: 

B = L (0) -|X|, 

S = L (0) . (3) 

It follows from (3) that B = L °' — |Xo| — J" sgn(X) d X and S t = |Xo| + sup{B s : s ^ t}, so the running supremum 
does not start at Bo if |Xo| > 0. 

Evidently it suffices to exhibit successful coupling strategies for (B, S); off the line X = 0, this Levy transform 
forms a 2 : 1 representation of the original pair (X, L' '); the two pre-images under the Levy transform meet 
together when the Brownian motion X hits 0. 

2.2 The reflection / synchronized coupling 

for immersed coupling of Brownian motion together with local time 

The above considerations show that it suffices to exhibit a successful immersed coupling between (a) the 
pair (B,S) of Brownian motion B and its running supremum S and (b) a copy (B,S) started with different 
initial conditions. Were the corresponding X and X not to agree at coupling, one could simply continue with 
synchronized coupling until |X| = |X| were to hit 0. However the reflection / synchronized coupling given below 
actually terminates with B = S and B = S, so at the end of this coupling we already have |X| = |X| = 0. 
Without loss of generality, suppose that Bq = L °' — |Xq| ^ Bq = L °' — |Xo|. 
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Definition 3 (Reflection / synchronized coupling). The reflection / synchronized coupling algorithm con- 
sists of two stages: 

1. Reflection coupling (dB = — dB) till the time Ti = inf{t : B t = B t } (the first time that B and B meet); 
then (if (B, S) is not already coupled with (B, S)) 

2. Synchronized coupling (dB = +dB), run from time Ti, first till the later time T2 = inf{t > Ti : B t = 
B t = So V So} that B = B first hits the level So V So after time Ti, then till the time T 3 = inf{t > T 2 : 
Bt = B t = Sj 1 V So} that B = B first hits the higher level Sj 1 V So after time Ti. 

Note that at the end of stage 2 we have B = S and B = S, so |X| = |X| = 0. 

It is possible for the coupling of (B, S) and (B, S) to be abbreviated to a one-stage (reflection) coupling in 
case So = So, for if it happens that B (and therefore B) both stay below So = So up to time Ti then coupling 
will be successfully achieved at time T COU pie = fi < T3. However we will see below that this case can be 
viewed as singular, as a consequence of Lemma 7. Moreover if Xo and Xo are of opposite sign then they will 
not couple at this stage: it will be necessary to proceed to completion of the synchronization stage so that 
By 3 = Sj 3 = Bj 3 = Sj 3 and therefore Xy 3 = Xj 3 = 0. 

If not abbreviated at the end of the reflection stage then the coupling will succeed at the time T coup i e = T 3 ; 
at that moment in time it is the case that simultaneously B = B (since they are coupled by synchronization 
after meeting at time Ti) and S = S(= B = B). The synchronized stage has been divided into two steps 
(over time intervals [Ti , T2 ) and [T2 , T3] ) since the [T2 , T3] step depends on the behaviour of B over the initial 
time interval [0,Ti]. Indeed, note that by construction (and particularly by choice of initial conditions Bo = 
L °^ — |Xo| ^ Bo = L °' — |Xo|) it is the case that B will stay below or equal to B until time T3, and hence 
Stj V So = Sy 1 V Sj 1 . The construction is illustrated in Figure 1. 
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Figure 1: Illustration of a successful reflection / synchronized coupling of Brownian motion B together with its 


running supremum S. 



The coupling is almost surely successful, since both the first and second stages correspond to times taken for 
real Brownian motion to hit specified levels. Indeed, the coupling of (B, S) and a copy (B, S) is equi-filtration, 
not just immersed, because the stopping times Ti, T2, T3 can be re-written as hitting times for B in its natural 
filtration (successively, from Bo to |(Bo + Bo), then from |(Bo + Bo) to So VSo, then from So VSo to Sjj V So), 
and similarly also as hitting times for B in its own natural filtration. (In particular, Ti can be rewritten as 
the hitting time of B moving from Bo to |(Bo + Bo).) The corresponding immersed coupling of (X, L' ') with 
(X, L' ') cannot be an equi-filtration coupling, because the natural filtration of X has to be augmented in order 
to supply appropriate randomness for the signs of the excursions of X from zero. (See Emery, 2009, Lemma 5 
for a similar augmentation in the more complicated case of BKR diffusions.) 
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There is a natural reformulation of reflection / synchronized coupling in terms of stochastic calculus: set 
d B = J d B up to the coupling time T2 , where the predictable control J is given very simply by 



Jt 



-1 for t < Ti (reflection stage) , 

-1 for Ti ^ t ^ T3 (synchronized stage) . 



(4) 



The failure of mutual immersion for the coupling of (X, lJ ') with (X, lJ ') is immediately apparent from the 
relevant stochastic differential equation 



dX 



sgn(X)Jsgn(X)dX, 



(5) 



which is an instance of Tanaka's classic example of a Brownian motion X, defined as a weak but not strong 
solution of a stochastic differential equation driven by a second Brownian motion J" Jsgn(X) dX. 

In the next subsection we will discuss optimality of the reflection / synchronized coupling amongst all 
immersed couplings. It is therefore helpful for purposes of comparison to supply a second successful immersed 
coupling, and so we conclude this subsection by outlining a different coupling which coordinates successful 
couplings of B with B and of S with S so that B does not couple with B before S couples with S, except of 
course for the case B = B . 

Definition 4 (Coordinated coupling). Suppose without loss of generality that Bo < Bo < (So V So)- The 
coordinated coupling is formed from an infinite sequence of consecutive cycles, each consisting of two phases: 

1. A phase of reflection coupling, till the following ratio is reduced by a factor of |: 



B-§' 
(S VS) - 



B + B) 



B-B 
2(SVS)-B 



(6) 



2. Then a phase of synchronized coupling till the ratio is increased by a factor of 2. 



The coordinated coupling is illustrated in Figure 2. Note that the coupling construction ensures that B > B 
over the sequence of cycles. Elementary estimates based on Brownian hitting times show that each individual 
cycle is almost surely of finite duration. 




Figure 2: Illustration of successful coordinated coupling of Brownian motion B and its running supremum S. The 
longer vertical bars on the time axis indicate switches between reflection and synchronized couplings. 



The ratio remains positive, and therefore B < B, throughout the sequence of cycles. Moreover, the initial 
ratio at the start of the n th cycle is 2~ n times the initial ratio at the start of the first cycle. The ratio remaining 
positive, it follows that the ratio at the start of the n th cycle converges to zero as in 00. Consequently 
liminf(B — B) — >• over this infinite sequence of cycles. 
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On the other hand |(B + B) remains constant in the reflection phase, while |(B — B) remains constant in 
the synchronized phase. Prom this we may deduce that (SVS) — |(B + B) halves at each successive termination 
of a synchronized phase. Outside synchronized phases, this quantity is non-decreasing. It may increase in the 
reflection phase, but only if B hits SVS, and this happens with probability given by the initial value of the 
ratio (6) at the start of the cycle, independently of previous cycles. The initial ratio decreases geometrically 
over successive cycles, so the first Borel-Cantelli lemma shows that the number of cycles in which B hits SVS 
is almost surely finite. Accordingly the sequence of values (S V S) — |(B + B) (considered sequentially at the 
start of each cycle) must almost surely decrease geometrically to zero after a finite initial number of cycles has 
taken place. Since B < B < (S V S), it follows that liminf ((S V S) - B) = and liminf ((S V S) — B) = when 
these liminf s are taken over the sequence of cycles. 

Finally, the above Borel-Cantelli argument can be adapted to show that there is a positive probability of 
the sequence of cycles being completed without SVS increasing from its initial level So V So- In this case the 
sequence of cycles must complete within the (almost surely finite) time interval during which B travels from 
Bo to So V So- Thus the total duration of the infinite sequence of cycles is finite with positive probability, and 
so we can use sample-path continuity to replace liminf by lim. Application of the Kolmogorov : 1 law shows 
that the total duration must therefore be almost surely finite, and hence that the coordinated coupling almost 
surely succeeds in finite time. 

In particular, note that if Bo < Bo < So = So then there is positive probability of coordinated coupling 
succeeding by means of B and B coupling exactly when they simultaneously first hit the level So = Sq. 



2.3 Strict optimality amongst immersed couplings 

The reflection / synchronized coupling strategy is faster than all other immersed couplings, in the sense that 
it strictly minimizes 

P [T couple > t] 

simultaneously for all t > 0, except perhaps for the singular case of So = So (this singular case is discussed 
around the statement of Lemma 7 below). Equivalently the distribution of the coupling time T coup i e for 
any immersed coupling exhibits stochastic domination over the distribution of T COU pie for the reflection / 
synchronized coupling (except perhaps in singular cases); moreover this domination is actually strict for any 
other immersed coupling. 

Before stating and proving a theorem which asserts this optimality, we first establish some preparatory 
lemmas. The first one concerns the coupling of two Brownian motions on [0, oo) that are stopped when the 
first one of them hits 0. 

Lemma 5. Suppose the planar process (U, V) is composed of two Brownian motions which are related 
by an immersed coupling, and suppose that (U, V) is started at a point (Uo,Vo) in the interior of the 
quadrant {(u,v) : u ^ 0,v ^ 0}. Let T be the first time that (U, V) hits the boundary of the quadrant. 
Suppose it is desired to construct the coupling so that P[(Ut,Vt) = (0,0)] = 1. This is possible if and 
only if Uo = Vo and the coupling is the synchronized coupling. 

Proof. Using the formalism of Ito (1975) (see also Ikeda and Watanabe, 1981, Ch. III. 1, and particularly 
the development in Kendall, 2001) and the representation of immersed Brownian couplings given in (1), the 
general law of (U, V) under an immersed coupling produces dU 2 = dV 2 = dt, Drift dli = Drift dV = 0, 
and dlldV = Jdt for an arbitrary adapted integrand J e [—1,1] which can be viewed as the control for the 
stochastic control problem of maximizing the objective function P [(Uj, Vr) = (0, 0)]. 

Without loss of generality we may suppose that Uo ^ Vo > 0. Note that under reflection coupling (J = — 1) 
the probability of (U, V) hitting the diagonal {(u, v) : u = v} before time T is given by 

V 

<D(Uo,V ) = -. . 

|(Uo+V ) 

We extend the definition of O to the case V ^ U > by setting O(U , V ) = <D(V , U ), so that 



0(11, V) 



U V 



|(U + V)'|(U + V) 
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An application of Ito calculus shows that if U > V > then, under a general control J e [— 1, 1], 

Drift d«D(U,V) = - (U-V)(l + J)dt, 

and this is non-positive, and vanishes only when J = — 1. A similar result holds for V > U > 0. On the other 
hand, if U = V > and J = +1 then (U, V) stays on the diagonal, so that <5(U, V) then remains constant. An 
argument using the Ito-Tanaka formula for semimartingales thus shows that 0(U, V) is a supermartingale for 
all immersed couplings of U and V, and becomes a martingale only under the strategy "use reflection coupling 
till (U, V) hits the diagonal or the boundary, then use synchronized coupling till (U, V) hits the boundary". 
It follows that O(llo, Vb) is the maximum of P [(Uy, Vr) = (0, 0)] over all immersed couplings, and is attained 
only by using this strategy. The lemma follows. □ 

As a consequence of the Lemma, we can prove the optimality of the reflection / synchronized coupling in 
the special case when B = B . This allows us to restrict attention to immersed couplings which preserve the 
ordering of B and B. 

Lemma 6. Consider the reflection / synchronized coupling (Example 3) for the special case B = B . 
This is the only optimal coupling amongst all immersed couplings started with B = B , and so (since 
the reflection stage succeeds immediately) the uniquely optimal way to proceed is to cease immediately 
if So = So, and otherwise to conduct a synchronized coupling of B and B until B = B hits So V So- 

Proof. If So = So then coupling succeeds immediately and there is nothing to prove. Suppose without loss of 
generality that So > So- Coupling cannot succeed earlier than the first time T at which B hits So, and if we 
employ synchronized coupling then coupling will succeed at this hitting time. 

This shows that synchronized coupling is optimal, but we require strict optimality. Consider a second 
coupling which does not employ synchronized coupling throughout. It then follows that there must be a 
moment, before time T, at which either B > B or B > B. We can apply Lemma 5 to U = So — B and 
V = So — B; it follows that if synchronized coupling is not employed right up to time T, then there is a positive 
probability that one of two possible cases has occurred: either B has already hit So by time T, or B has not 
yet hit So by time T. In the first case, the properties of Brownian motion B show that almost surely Sy > So, 
and so successful coupling must occur after B travels from So to S^. In the second case, successful coupling 
must wait at least until B and B meet after time T. 

It follows that, for any coupling other than synchronized coupling, (a) the coupling time can be no less than 
T, (b) there is a positive chance of it being strictly greater than T. This establishes the required stochastic 
domination (since T is the hitting time of a real Brownian motion started at B = B and rising to So > So VSo) 
and so the lemma follows. □ 

In passing, we are now able to explain the reason why it is appropriate to describe as singular the case 
when Bo, Bo < So = So- In this case it is possible for full coupling of (B, S) with (B, S) to succeed as soon as B 
first meets B, so long as B and B do not hit So = So (as noted in Section 2.2, this need not imply success of 
the coupling of X with X at that time) . The next lemma shows that if So 7^ So then this early success cannot 
occur. 

Lemma 7. Suppose So 7^ So- Then an optimal immersed coupling of (B,S) and (B,S) succeeds exactly 
at the first time when B, S, B, and S simultaneously coincide. 

Proof. Consider first the case Bo = Bo- As shown by Lemma 6, it is then the case that the only optimal 
immersed coupling is provided by synchronized coupling until B — B first hits So VSo, and the characterization 
of coupling by simultaneous coincidence is immediate. 

Consider the case Bo < Bo (the case of Bo < Bo is entirely similar). It is a consequence of Lemma 6 
that optimality of the immersed coupling implies that the relationship B ^ B must persist till full coupling is 
successful. 

So further suppose that So < So. In that sub-case it follows from B ^ B that the relationship S ^ S must 
persist till full coupling is successful. Coupling cannot succeed till S hits S, and when that happens we must 
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have B = S. But in this sub-case we also have B ^ S ^ S and B ^ B ^ S. Consequently full coupling must 
succeed when S first hits S, and at that time B, S, B, and S simultaneously coincide. 

Suppose on the other hand that So < So- In that sub-case again, full coupling cannot succeed before S hits 
S, at which time it is necessary that B also hits S. If it is further the case that B = B at that time, then full 
coupling succeeds in the manner prescribed by the lemma. If on the other hand B < B at that time then (by 
the properties of Brownian motion) there are instants immediately after this time at which B < S < B < S, 
and we can proceed as above. □ 

We shall now show that the distribution of the coupling time T coup i e under any immersed coupling can 
be dominated in the limit (as N — > oo) by the distribution of T COU pie under an immersed coupling whose 
predictable control J is restricted to values ±1 (thus, a "bang-bang" control), and moreover such that J is 
constant on stochastic intervals [t^ ,tL \) defined as follows. For any positive even integer N > 0, consider 
the one- dimensional lattice £' N ' (depending implicitly on Bo, and Bo) 

£ (N) = B + B ° N B ° Z = {Bo + ^(B -B ) : k = 0, ±1, ±2, . . .}. 

We define a mesh, a sequence of stopping times = t N ^ < t[ N ' < Tj N ' < . . ., as a sequence of "crossing 
times" for this lattice: 

tjK = inf{t>x« N) : B t e£( N '\{B T « N) }}. 

Sampling using this mesh of stopping times has the effect of discretizing the Brownian motion B into a random 
walk with steps ±^(B — B ). 

Note, for an immersed coupling restricted to a control J which is locally constant with J = ±1 on each 
stochastic interval [t^ , of the mesh: 

1. both B and B, when sampled at times = x Q N ' < t[ N ^ < < . . ., belong to the lattice £' N ', since B 
is obtained from B using a predictable control J formed from synchronizations and reflections and which 
alters only when B belongs to the lattice; 

2. because N is even, T coup i e belongs to the set {t N ', t[ , Tj , . . .}; 

3. finally, our candidate for optimality, the reflection / synchronized coupling (Example 3), can itself be 
viewed as one of these couplings, since the control changes from +1 to —1 exactly at one of the stopping 
times in the mesh. (This is the reason why it is convenient to work with meshes of stopping times, rather 
than decompositions of the time axis into disjoint dyadic intervals.) 

We can now summarize and prove a result stating that an optimal immersed coupling can be approximated 
in distribution by appropriately chosen "bang-bang" controls of the above form. The proof is related to the 
method of proof of Emery (2005, Proposition 2); however here we need the control J to have the "bang-bang" 
property rather than simply to be locally constant, and to be composed of stopping times drawn from a mesh 
of stopping times as specified above. 

Lemma 8. For any fixed t > 0, any optimal immersed coupling of (B,S) and (B,S) can be approxim- 
ated weakly over [0,t] (when viewed as a probability distribution on the metric space of A- dimensional 
continuous trajectories, equipped with the sup-norm) by "bang-bang" immersed couplings for which the 
control J takes values ±1 only, and only changes at hitting times belonging to some mesh. 

Note that the lemma does not assert that the "bang-bang" couplings are successful! 

Proof. Consider a general immersed coupling determined byB = Bo + /JdB and subject to the constraint 
that the coupling is synchronized once B and B have met. (By Lemma 6 all optimal immersed couplings must 
be of this form.) Since |J| ^ 1, for each t > we have E J* J 2 ds < oo, and moreover for each e > we 

may find continuous predictable f with E J |f— J| 2 ds < £ 2 /4 (for example, f t = -p- J*_ 6 (t — s)J s d s for 
sufficiently small 5). It then follows that for sufficiently large N we may approximate f in L 2 by piece-wise 
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constant j' c ' such that j' c ^ e [— 1, 1] is predictably constant on each dyadic interval [tJ. N ' , ) of the mesh 



and E 



j^|f-jM| 2 ds 



< £ 2 /4, hence 



E 



I J — J 



[c]|2 



ds 



< £ 



Doob's submartingale inequality then implies that we can control 



sup 

s<t 



JdB 



J lc| dB 



so that Bo + J j' c ' d B is a good path- wise approximation to B = Bo + / J d B. 

While j' c ' is piecewise-constant on stochastic intervals related to the mesh, it does not take values in {±1}. 
We need an approximation based on a "bang-bang" control j' bb ', which is constrained by J^ bb ' G {±1} as well 
as by the requirement that j' bb ' is predictably constant on stochastic intervals [t{, M ', ) which now must 
be formed on a new mesh, defined for some still larger even integer M = 2 r N, for an integer r > 0. Given 
M > N, we define jI bb )i( M ) to "track" J [cl in the following co-adapted way: 



r [bb];(M) 

I (M) 



-l if Jo* j[ bb ];( M >du< f k j^du, 

-1 if J^ M> j[^];(M) du> j^ M) jlc] du . 



Since |J| ^ 1, it follows that we have the following bound for s € [0,t]: 



J [cl du- 



[bb];(M) 



du 



2 sup {(t{™ At) 



JM) 



At) : k = 1,2, 



(7) 



(8) 



converging almost surely to zero as 2 r = M/N — s- oo. 

Consider the sequence of two-dimensional processes {(B,Bo + J" j' bb ^ M ' d B) : M = 2 T N}, defined on 
the time-range [0, t]. The one-dimensional coordinate processes being Brownian motions, it follows that this 
sequence is tight. Any convergent subsequence converges to a limit for which the one-dimensional coordinate 
processes are Brownian motions; moreover, using (8), we may deduce that in the limit the product of the 
pair of one-dimensional coordinate processes is equal to the sum of a martingale and the integral J^ c ' du. 
Hence by semimartingale Ito calculus the limit has the law of (B, B + J j' c ' d B), no matter what convergent 
subsequence is chosen, and therefore by the theory of weak convergence we may deduce that the sequence of 
random paths (B, B + J J [bb)i(M) d B) converges weakly to this limit. 

It follows that we can choose a sequence of "bang-bang" controls j' n ', constant on appropriate meshes 



(M„ 



l k+l 



k = 1,2,...} (with M n — > oo), such that (B,Bq + Jj'^'dB) converges weakly (using 



supremum norm over the time interval [0, t]) to the immersed coupling (B,B) which was originally under 
consideration. □ 

We can now argue for optimality of the reflection / synchronized coupling in the general non-singular 
case (So ^ So). We need only consider the case when B 7^ B , since Lemma 6 covers the case of B = B . 
From Lemma 6 it suffices to consider only those immersed couplings which are constrained to be synchronized 
couplings once B and B have met. Employing the terminology of subsection 2.2, we set Tj = inf{t : B t = B t }. 
Thus we need consider only those immersed couplings for which J t = 1 once t > Tx- 

Theorem 9. Suppose that Bo ^ Bo and So ^ So- The reflection / synchronized coupling (Example 3) is 
optimal amongst all immersed couplings of Brownian motion together with local time. 

Proof. Note that by Lemma 6 we may restrict attention to immersed couplings for which (without loss of 
generality) B ^ B, and such that B = B after Ti = inf{s : B s = B s }. Moreover, by the argument of Lemma 7, 
at the coupling time T couple we must have B Tcouple = S T „ upIe = B Tcouple = S Tco „ ple . 

The first step is to use the weak approximations (B, Bo + J j' n ' d B) of (B, B) (as given in Lemma 8) to build 
successful immersed couplings of (B,S) and (B,S) with coupling times which are in the limit stochastically 
dominated by the coupling time derived from (B, B). For convenience we employ the Skorokhod representation 
of weak convergence; augmenting the probability space if necessary, we construct a copy (B' n ',B' n ' = Bo + 
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J J*'<^ dB in) J of (B, B + J J (T1) d B) on the same probability space as (B, B) such that almost surely B (n) -> B 
and B' n ' — > B uniformly on the time interval [0, t]. (We note in passing that this construction need not respect 
the underlying filtration. The stochastic integrand J*'' n ' and the stochastic integral J" J*>( n ) dB' n ' are defined 
with respect to the natural filtration of B' n ', which need not immerse in the original filtration!) 

Although the target coupling of (B, S) and (B, S) couples at T coup i e , we should not suppose that (B' n ', S' n ') 
and (Bf n ',S' n ') couple at this time (using S' n ' and S' n ' to denote the corresponding supremum processes). 
However we can modify (B* n ', B' n ') to produce a coupling which does not succeed much later than the original 
coupling. 

Indeed, we have restricted attention to immersed couplings such that at the coupling time T coup i e we have 
B T couple = S Tcouple = B Tcouple = S TcoupIe . Accordingly we may choose a sequence £ n -> such that 



3' n) , S[ n) , B' n) , S' n) all lie within ±e n of each other 

I couple ' couple ' couple ' couple 



> 1 



Accordingly, if we set 

T 3 tn) = inf js : B< n) , S< n) , B< n) , S l s n) all lie within ±e n of each other} 



then 



> > t 



«S P [T couple > t] + En . 

But at time T^ n) we can modify the construction of (B (n) ,S (rt) ) and B (n) ,S (n) ) to use the reflection / 
synchronized coupling (Example 3), obtaining successful coupling at time T^ ple > Tg . As £ — > so we can 

deduce by scaling that the extra time T c ^ le — Tg required for success of this final coupling must tend to 
zero in probability. 

It follows from these arguments that the infimum of the probability of failing to couple before time t, for 
any fixed t > 0, 

P [Tcoupie > t] , 

can be approached by considering P [T COU pie > t + £rJ for suitable £ n — > and immersed couplings based on 
"bang-bang" controls j' bb ' constrained by change only at stopping times taken from meshes = t N ' < t[ < 
Tj N ' < . . ., and which become synchronous after B and B first meet. 

Consider such an immersed coupling with control j' bb ^. For a fixed t > 0, we consider the following value 
function, defined for ^ u < t: 

V(u;b,b, s,s) = 

P T C o U pie > t — u | B u = b, B u = b, S u = s, S u = ?; reflection / synchronized coupling . (9) 

We are particularly interested in the discrete-time process obtained by sampling at stopping times taken from 
the specified mesh, but stopping at the terminal time t: 

{z u =V(^>At|B T , N]At ,B T , N)At ,S T « N)At ,S T , NlAt ) : n = 0,l,2,...}. 

It follows by definition that V(u; B u , B^, S u , S u ) is a bounded martingale under the reflection / synchronized 
coupling, and hence under this coupling Z is a discrete-time martingale since it is obtained from the bounded 
process {V U At : u ^ 0} by sampling at stopping times. We shall show that Z is a supermartingale under the 
coupling specified by J^ bb ^, and moreover that the martingale property cannot hold if J^ bb ^ = +1 over the 
initial time interval [0,t^ N '). Arguing inductively, this suffices to establish the theorem. 

The crux of the matter is to consider the behaviour of the value function at time zero if the initial segment 
of coupling is synchronized. To this end, we make a special construction of the reflection / synchronized 
coupling referred to by the value function: we suppose two independent Brownian motions are employed (both 
begun at 0), namely B' r ' to drive the reflection stage of the coupling, and B' s ' to drive the synchronized stage. 
We set 

. xi N;s) =inf{t>0:|B[ s) | = i(B o -B o )}, 
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T w>t 



• j[ to be the time when B^ first hits — | (Bo — Bo), corresponding to the end of the reflection stage, 

( \ (r) fr) 

• and M lr) = sup{B^' ' : s < Tj' '} + B to be the maximum level achieved during the reflection stage. 

Then the reflection / synchronized coupling time corresponds in law to T3*' = Tf +inf{s : E>!; S ' + |(Bo + Bo) = 
max{S V So, M (r) }, and so Z = 

So consider the effect of commencing with a session of synchronized coupling. We consider two possible 
cases. Suppose in the first case that 

M (r) -^(B -B ) < S0VS0. 
Then we can represent the initial session of synchronized coupling by using B* s '|, (n; S ),, and then replacing 

(si fsl (si (*1 

B| by B (N j — B (N|S) . Evidently the distribution of T3 is unaffected by this change. If furthermore 

t + Tj ' T 1 

M M + ^(B -B ) s$ S0VS0 

then the reflection stage will start at time t[ N,s ' at level b' s (N;s) , moreover by the end of the reflection stage 

the supremum of the coupled processes will not exceed So V So, and the subsequent synchronization stage will 
have to move from B ( ^ N . S) + |(B + B ) to So V So- It follows that T3* 1 is still the coupling time. If on the 
other hand 

M (r) + ^(B -B ) > SqVSo, 
then there is a possibility that the supremum of the coupled processes will exceed So V So- However the 

(si 1 — — 

subsequent synchronization stage will still have to move from B IN sl + ^(Bo + Bo) to So V So, but may have 
to move even further. Thus the coupling time still cannot occur earlier than T 3 . 



Suppose in the second case that 



M (r) -l(B -Bo) > SqVSo. 



Then we can represent the initial session of synchronized coupling by using an independent copy BL .(n |S ), of 



[0,T 



and restarting the construction at the new starting points Bq ± ^ (Bo — Bo), Bq ± ^(Bq — Bq) (same 



sign for each initial increment). Regardless of the sign of the initial increment, coupling occurs 

f* I 

so is delayed relative to T3 . 

It follows from these arguments that an initial session of synchronized coupling followed by reflection / 
synchronized coupling cannot increase the probability of successful coupling by time t compared with that of 
reflection / synchronized coupling, and has a positive chance of reducing it; consequently Z is a supermartingale 
and the martingale property for Z cannot hold if J^ bb ^ = +1 over the initial time interval [0,t[ N '). 

Since Z is a martingale for the reflection / synchronized coupling, this suffices to establish the theorem. 

□ 

In the singular case the value function (9) takes on a more complicated form, and it is no longer clear 
whether or not the reflection / synchronized coupling is optimal amongst immersed couplings. 



2.4 Rate of optimal immersed coupling 

This subsection elicits the rate at which the reflection / synchronized coupling occurs. This is accomplished 
by calculating the moment generating function of the optimal immersed coupling time in the non-singular case 
of Sq So, supposing (without loss of generality) that B > B . 
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Theorem 10. Let J CO upie be the coupling time for Brownian motion together with local time, equivalently 
for Brownian motion B together with its swpremum S, using the reflection / synchronized coupling 
described above. Let B and S be the corresponding coupled quantities. Under the non-singular conditions 
So = Bq > Bq = So, the coupling time has the following moment generating function: 



E [exp(— oCY coup ie)] = 1 -sinh ( B o - B )^ logcoth 




Proof. Recall the notation of subsection 2.2, and bear in mind the stipulation that So 7^ So. It is required to 
calculate the moment generating function 

E [exp(-cxT COU pie)] = E [exp(— <xT 3 )] . 

We can express T 3 as the sum of (a) the Brownian hitting time Ti, being the time taken for B to pass from B 
to I (B + B ) , (b) a further Brownian hitting time 1 2 — Tl, being the time taken for B to pass from | (B + B ) 
to So VSo, and finally (c) a randomized Brownian hitting time T3 — T2, being the time taken for B to pass from 
So V So to Mi = sup{B t : t ^ Ti}. Amalgamating the second two of these, 

T 3 = Ti+H 1 fmax{S VSo,Mi}-i(Bo + B, 

where H 1 (a) is the time taken for a standard Brownian motion to pass from to a. 

We outline the calculations for the special case Bo = So and Bo = So (though the calculations can be 
extended to the general case). Thus we are concerned with the moment generating function of Ti + H 1 (Mi — 
|(B + B )). 

The first task is to investigate aspects of the joint distribution of Ti and Mi, specifically 

Q a (a) = E[exp(-ocTi) ; Mi < a + B ] . 

This can be calculated using excursion theory, for example by adapting the calculations of Rogers and Williams 
(1987, §56). For convenience we set a* = v2a and b = Bo — |(Bo + Bo) = §(Bo — Bo). Suppose that excursions 
are marked using an independent Poisson(a) point process on the time axis. Then we distinguish the following 
kinds of excursions of B from Bo, noting the rates at which they happen when the excursions are viewed as 
points of a Poisson process of excursions with respect to local time: 

1. upward excursions which do not rise above the level q+Bo but are marked (occurring at rate |(a* coth(<x*a) 

2. upward excursions which rise above the level a + Bo (occurring at rate ~); 

3. downward excursions which do not fall below the level |(Bo + Bo) but are marked (occurring at rate 

i(a*coth(a*b)-i)); 

4. downward excursions which fall below the level §(Bo + Bo) (occurring at rate ^t); 

5. downward excursions which fall below the level | (Bo + Bo ) but are not marked before hitting | (Bo + Bo ) 
(occurring at rate ^- cosech(a*b)). 

These rates are computed as in the discussion of Rogers and Williams (1987, §56), based on the identification 
of the law of the Brownian excursion discussed there. Thus Q«(a) can be computed as the probability that 
we see a downward excursion falling to level |(Bo + Bo) before it has been marked (that is, an excursion of 
type 5) before we ever see upward excursions rising to level a + Bo (of type 2), or staying below this level but 
marked (type 1), or downward excursions which do not fall below the level |(Bo + Bo) but are marked (type 
3), or downward excursions which fall below the level |(B + B ) but are marked before hitting |(Bq + B ) 
(type 4 but not type 5). 
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We can therefore use Poisson point process theory to compute 

Qa(a) = 



^- cosech(a*b) 



i + §(a* coth(a*a) - i) + §(a* coth(a*b) - £) + (i) 

cosech(a*b) 
coth(a*a) + coth(a*b) 



sinh(a*a) 
sinh(a*(a + b)) 



(10) 



Consequently the desired moment generating function is given by 



E 



exp(-a(T 1 + H 1 (M 1 --(Bo + B ))) 









E 


exp - 


Jo 
















E 










Q«(da) 



a* sinh(a*b) 



'°° E [exp (-aH x (a + b))] 
o sinh 2 (oc*(a + b)) 



da 



(once again using b = |(B — B )). We use the reflection principle 



[H 1 (a + b) >t] 



□ +b 



- 2 / 2 d1 



and integration to deduce the moment generating function for the optimal immersed coupling time when 
So = Bq 7^ Bo = So, 



E [exp (-aT COU pie)] 



E 



exp ( -«(Ti + H 1 (Mi - -(Bo + B ))) 



a' 



l-sinh( — (B -B ))logcoth — (B -B 



a 



(11) 

□ 



These excursion-theoretic calculations have been checked by simulation, although direct simulation is com- 
putationally demanding because of the heavy tails of the Brownian hitting times involved in the reflection / 
synchronized coupling. This can to some extent be mitigated by the use of "Rao-Blackwellization" using the 
known formula for the moment generating function of the Brownian first passage time. 



2.5 Comparison with maximal coupling 

It is natural to ask whether the optimal immersed coupling is in fact a maximal coupling. Numerical evidence 
is that this is not the case, as is readily seen by computing the total mass of the minimum of the joint densities 
for (B t ,S t ) and (B t ,St), and comparing with the coupling probability for optimal immersive coupling. Prom 
Revuz and Yor (1991, Ex 3.14 part 2°) the joint density of B t and S t (supposing So = Bo = 0) is given by 

2~2s-b / (2s-b) 2 \ „ , , 

exp — forbsCs,s>0. (12) 



7rt t V 2t 



(This follows readily from the reflection principle.) We can use this to evaluate numerically the moment 
generating function of the maximal coupling times (based on the maximal coupling derived from the method 
suggested by Sverchkov and Smirnov, 1990). Indeed, if T max i ma i is the coupling time for the maximal coupling, 
with So = Bo = and So = Bq = —bo < 0, then 



I 5 [Tmaximai ^ ti 

' 2 



2s -b (2s-bi2 . f 2s-b+b / bg b (2s-l, , . , , , 
-^e « nun U — - exp f 1 )^dbds. 
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Changing variables from b ^ s, s ^ = max{0, —bo} to ^ v = 2s — b< oo, 0^u = s^v = 2s — b, and 
integrating out u over the range (0,v), we obtain 



P[T, 



maximal 



<t] = 



— . , . v + b / bg + 2b v > , ( , 

T e ='mmjl, — exp — | )■ dv, 



and thus we obtain a double integral for the moment generating function of the maximal coupling time: 

E [exp(-aT maxima i)] 



2 r°° av 2 e-' xt -^ 



min < 1 



v + b 



exp 



2b v 



2t 



dvdt (13) 



A numerical comparison of the moment generating functions for the two couplings is given in Figure 3: it is 
evident that the two do not agree. This is numerical eviedence that the optimal immersed coupling is distinct 
from (and therefore slower than) the maximal coupling. 




Figure 3: Comparison of moment generating functions for coupling times of (a) optimal immersed coupling (dashed 
curve) and (b) maximal coupling (continuous curve) of Brownian motion B and its running supremum S. Initial 
conditions are B = S = and B = S = — 1. 



3 Mutually immersed couplings of Brownian motion together with 
local time 

As noted in subsection 2.2, the reflection / synchronized coupling is an equi-filtration coupling when viewed 
as a coupling for the absolute value |X| of Brownian motion and L' ' its local time at 0. However it is 
not equi-filtration when the absolute value |X| is replaced by X, the Brownian motion itself. This follows 
immediately from consideration of the stochastic differential equation (5), which is of Tanaka type when 
viewed as generating the coupled Brownian motion X from Jsgn(X)dX. Unless X are identical, it is not 
possible to extract statistically appropriate signs for the excursions of X from the natural filtration {!F t : t ^ 0} 
of X. 

The coupling of (X, L ,0) ) with (X, L (0) ) can be modified to be equi-filtration by replacing sgn(X t ) and sgn(X t ) 
in the coupling control by time-delayed versions, at the price of delaying the time of successful coupling. Given 
a positive non-increasing function 4> : [0,oo) — > [0,oo), we introduce the delayed time-change 

o-(t) = t-(i|)(t) At). (14) 

For t > we define a new coupled Brownian motion X t , starting at Xo = Xo but defined up to time T3 as the 
solution to a time-delayed version of (5): 

dX t = sgn(X (T ( t) )J t sgn(X ff ( t) )dX t . (15) 
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Here we use (4) to define the control J in terms of X via the stopping time Ti : 



= j-1 ift^Tj 
1+1 otherwise. 

Of course this exploits the remark in subsection 2.2, that Ti and T3 can be defined as hitting times for 
B = L' ' — |X|. Adopting the usual convention that sgn(O) = 1, and arguing from the positivity and continuity 
of ij), we can solve the stochastic differential equation (15) step by step over successive small time intervals. This 
ensures that sgn(X CT ( t )) and sgn(X cr ( t )) are defined and measurable with respect to S^ft) = o~{X s : s ^ <r(t)}. 

This construction makes it clear that the coupled X is immersed in the natural filtration of X. However the 
reverse also holds: 

Lemma 11. For X defined in terms of X using the time- delayed stochastic differential equation (15), 
it is the case that X is immersed in the natural filtration of X, so that the coupling of X and X is 
equi- filtration. 

Proof. The control J appearing in (15) satisfies J = — 1 up to the stopping time Ti. Accordingly X = X* up to 
Ti, where 

dXt = -sgn(X* a(t) )sgn(X ff(t) )dX t . 
However we may re-write this last stochastic differential equation as 

dX t = -sgn(X ff[t) )sgn(X* ff(t) )dXt, 

and so we find that 

dX t = -sgn(X (T (t))sgn(X (T (t))dXt 

holds up to time Ti. Since Ti is a hitting time of X, it follows that X stopped at time Ti is adapted to the 
filtration of X, and thus that Ti and thus J are adapted to the natural filtration of X. 
Arguing from time Ti onwards, since J t = 1 for t > Ti, we can re- write (15) as 

dX t = sgn(X (T(t) )J t sgn(X lT(t) )dXt = sgn(X (T(t) ) sgn(X a(t) ) d X t for Ti < t < T 3 . 

It follows that X up to time T3 is adapted to the natural filtration of X. This establishes the mutual immersion 
property. □ 

Of course X 7^ X, and therefore we cannot assume that the equi-filtration coupling will have succeeded by 
time T3. However we can use the properties of the reflection / synchronized coupling to argue not only that 
the paths of |X| and |X| will be close in probability, but also that the same is true of the respective local times 
L< 0) and L«». 

Lemma 12. There is a sequence of equi-filtration couplings of (X' n ', L'°' ; ' n ') with (X, L' ') such that 
(IX^'I, L'°' ; ( n ') converges in probability under supremum norm to (|X|,D ') (using the reflection / syn- 
chronized coupling): 



sup { |X[ 
t i- 



X t 



+ 



f (0);(n) _r(0) 



>4" 



< 



It suffices to exhibit a sequence for which the supremum converges to in probability. Note that we do not 
assert that X' n ' converges to X in probability: were this true, it would contradict the known fact that the 
stochastic differential equation (5) of Tanaka type can have no strong solutions. 

Proof. To simplify notation, we take Xo = 0. Note also that by construction t °' = L °' and Xo = Xo. 
First note that, if we set S = L (0) and B = L (0) - |X|, then 

dB = sgn(X)sgn(X ff )Jsgn(X (T )dX. 
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By choice of initial conditions, B = B . It suffices to show convergence in probability of sup t {|B t — B t |}, the 
supremum norm of the difference. Because |J| = |sgn(X)| = |sgn(X)| = 1 we can use Doob's L 2 submartingale 
inequality, the L 2 isometry for Brownian stochastic integrals, and Jensen's inequality to deduce that 



E 



sup{(B t -B t ) 2 } 
t 









•t 


= E 


supj 


( 



















oo 




4E 
















(sgn(X) sgn(X a ) J sgn(X a ) - J sgn(X)) d X 
| sgn(Xt) sgn(X (T(t) ) sgn{X a[t) ) - sgn(X t )| 2 d t 



E[|sgn(X t )sgn(X ff(t 
= 32 



1| 2 1 dt + 8 



Er| S gn(X CT(t ))-sgii(X t )| 2 ldt 



P [sgn(Xt) ^sgn(X ff(t) )] dt + 32 P [sgn(X ff(t) ) ± sgn(X t )] dt. 

Jo 

Both X and X are Brownian motions, though not necessarily begun at 0. It is therefore immediate that 
P [sgn(X t ) ^ sgn(X a(t) )] and P [sgn(X a(t) ) ^ sgn(X t )] are both dominated by P [sgn(X t ] ^ sgn(X ff(t) )] for X 
a standard Brownian motion begun at (since on average X will have to travel further to change sign than the 
standard Brownian motion X). Moreover rotational symmetry of the standard normal distribution reveals, if 
t>a|>(t), 



[sgn(X CT(t )) ^ sgn(X t 



= ^[PW^Xt-X^jl] = \ 



[sgn(X ff(t) ) + sgn(X a(t) + (X t - X a(t) ))] 

|X CTCt )l/Vt-il)(t) 



KXt-x^tjD/vARt; 



< 



' Mt) 
t-Mt) 



= - tan 

7T 



' Mt) 
t-Mt) 



Thus we obtain, for e > i(>(e), 



E 



sup{(B t - B t ) 2 } 
t 



s$ 64 



[sgn(X CT(t) ) ^sgn(Xt)] dt 

1 roo 



s$ 64 e 



71 



tan 



t-Mt) 



dt 



< 64 £ + 



71 



I Mt) 
t-Mt) 



dt 



For any e G (0, §), with 



Mt) = 



if t > £ , 



(t-£ + l) 3 

we find (a) \p is positive non-decreasing over [e, oo), (b) Mt) ^ t for t > e, and hence (c) 



E 



sup{(B t - B t ) 2 } 
t 



s$ 64 1 + 



71 



V(t/e)(t-e + l) 3 -e : 



dt £ 



< 64 1 + 



71 



/4 u+1 



du e = 105.557... x e . 



The result follows, because L 2 convergence of random variables implies convergence in probability. □ 

Theorem 13. There are successful equi-filtration couplings of Brownian motion together with local time 
starting from any pair of initial conditions (Xo,Lq ') and (Xo,Lq '). 

The proof makes it plain that the coupling time for reflection / synchronized coupling can be approximated 
arbitrarily well in distribution by the coupling times for suitable equi-filtration couplings. Consequently there 
can be no optimal equi-filtration coupling for Brownian motion together with its local time. 
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Proof. It suffices to show that there is a sequence e n — > such that the following holds for all reflection / 



synchronized couplings: if 



|X | - |X | 



(0) 



(01 



< £ n and X = Li ' then 



[Tcouple > 4- n ] S? 4" 



(16) 



For then we may construct a successful equi-filtration coupling as the concatenation of a series of equi- 
filtration couplings. In the first stage, Lemma 12 can be used to select an equi-filtration coupling approximating 
a reflection / synchronized coupling (continued up to but not including time T 3 , the end of the synchronized 



stage) such that we have the following control on the left-limits X, 



X t (d 

1 3 



f (0) 

-1 L T (1) 



1 3 



and L l 



1 3 



|X (ii I — |X (i) 

1 -P 1 P 



+ 



1 3 



— L 



(0) 

- T d)_ 

1 P 



> £i 



< 



Moreover, it follows from the construction of the reflection / synchronized coupling that |X m 

\x jW _\ = Z% = o. 

'3 lo- 
ll 



At time T, 



■ (0) 
- T CD_ 



X makes a small jump so that X T (i) 



X T (i) 



while X, L (0) , L (0) , X and L (0) trajectories 



remain continuous. The second and further stages are implemented by repeating this procedure. 

To be explicit, the construction continues through further stages n = 2, 3, . . ., such that at the end T ; 
of stage n, conditional on successful fulfilment of all previous stages, we have |X T cn) | = L 

I 3 



(0) 

- T ( n j 
1 p 



= |X, 



L tn] =0, moreover if n ^ 2 then T 3 



(n) 



. Ta (n-1) <4 _n 



and 



|X T tn)_| 

1 ■a 



|X T (n) | 
1 -a 



MOI 



"10) 



> £r 



< 



(17) 



(We suppress the conditioning on previous stages for the sake of simple notation.) Stage n is implemented (a) 
by using a reflection / synchronized coupling of (|X|, L' ') with (|X|, L* ') which succeeds before time 4~ n (this 
has probability 1— 4~ n ), and also (b) by invoking Lemma 12 to continue (X, L' ') by an equi-filtration coupling 
such that the maximum difference over this stage between the immersed coupling component (|X|, L' ') and the 
equi-filtration coupling component (|X|, L' ') is less than e n (with conditional probability at least 1 — 4~ n ). It 
follows that the conditional probability of the n th stage (n ^ 2) completing successfully is at least 1 — 2 x 4~ n . 
At the end of stage n we impose a small jump on X so that X, 



Thus, with total probability at least 



= X-pfn) 
1 3 



Y_ 2 x 4T 



n=2 



7 

12 



for all n, stage n is of duration less than 4 n , and at time T3 



(0) 



we have X_(n) = L 

T 3 T. 



= 0, and (X (n) , L 



0] 



can be approximated by the equi-filtration coupling component (X T (n),L |n) ) with maximum difference less 

' 3 T 3 

than e n . 

It follows that the bound (16) can be used together with Lemma 12 to construct an equi-filtration coupling 
which has probability at least ^ of succeeding in finite time. In case of default at any stage, one can then 
restart the sequence of couplings, so successful immersed coupling is almost sure to happen. 



It remains to establish the bound (16). Suppose that 



IXn 



|X 



(01 



"(0) 



< e and X = Li 



In 



case L< 0) VL< 0) 



fL, 



0) 



-0 I'M) 

generating function given by 



X |)V(L 



(0) 



|Xq|), we know from Theorem 10 that the coupling time has moment 



E[exp(-aT coup ie)] 



sinh 



b log coth 



a b 
2 2 



where b = |(Lq ' — |Xo|) — (Lq ' — |Xo|)| ^ £. Moreover by construction of the reflection / synchronized coupling 
the coupling time increases monotonically in b (note that it would be a faux pas to assert that monotonicity of 
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moment generating function implies stochastic monotonicity). Thus in this case the coupling time is dominated 
stochastically by a random variable with moment generating function 1 — sinh (\p§i) log coth i 
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In case L, 



(0) 



vl °> > 



ir-ixoi 



V (Lq — |Xq I ) , the reflection / synchronized coupling is completed by a 



passage of the synchronized Brownian motions from level (Lq ' — |X |) V (L °' — |X 



to L^VLj 01 , and this is 



dominated stochastically by the time taken for standard Brownian motion to move from to level e. 

The sum of both components will stochastically dominate any reflection / synchronized coupling time for 



which |X |-|X 



bound (16) follows. 



Ln 



'(0) 



< e and Xq = L ( 



10) 



Since both components tend to zero in distribution, the 

□ 



The following corollary will be of assistance when constructing equi-filtration couplings of BKR diffusions 
in the next section. 

Corollary 14. The equi-filtration coupling of Theorem 13 can be localized in the following sense: for 
fixed 5 > 0, for all sufficiently small e > and for all pairs of initial conditions (Xo,Lq ') and (Xo,Lq ') 
with 



|X | - IXol 



(0) 



Ln 



< e, 



we can construct a successful equi-filtration coupling of [X, L' ') and (X, L' ') such that 



one o/(X,L (0) ) and (X,L (0) ) does not stay within ball((X , L ), 6) 



< 



Proof. Prom the proof of Theorem 13, it follows that for sufficiently small e we can obtain uniformly arbitrarily 
small probability of the equi-filtration coupling failing to couple within an arbitrarily small period of time. 
The corollary then follows by observing that it follows from continuity of Brownian motion and Brownian 
local time that over a sufficiently small period of time we can ensure that the probability of there being large 
deviations either of motion or of local time is arbitrarily small. □ 



4 Application to BKR diffusions 

The BKR diffusion was introduced by Benes, Karatzas, and Rishel (1991) as part of an investigation into a 
control problem which possesses no strict-sense optimal law: it is a two-dimensional diffusion (X, Y) for which 
X and Y are Brownian motions connected by 

sgn(X)dX + sgn(Y)dY = 0. (18) 

Essentially (X, Y) diffuses linearly in a Brownian fashion along the boundary of a square {(x, y) : |x| + \y \ = £}, 
except that I increases according to a local time driving term whenever (X, Y) visits one of the vertices of 
the square. Emery (2009) studied filtration questions concerning this process, and in particular showed that 
the natural filtration of (X, Y) is Brownian even when (X, Y) is started from the origin. In this connection he 
asked a specific question (Emery, 2009, Question), which can be stated concisely as follows: given two initial 
points (Xo, Y ) and (X , Y ), neither of which is the origin, is there an almost surely successful equi-filtration 
coupling of BKR diffusions (X, Y) and (X, Y) started from these two points? An affirmative answer would lead 
to a constructive proof of Brownianity of the filtration of (X, Y) started from the origin. An immersed coupling 
is exhibited in Emery (2009, Lemma 5); however equi-filtration is required for the purposes of obtaining a 
constructive proof of the filtration result. 



4.1 Sketch construction of immersed coupling for the BKR diffusion 

The work of Sections 2, 3 suggests a direct strategy for constructing successful equi-filtration couplings of 
BKR diffusions; start with an almost-surely successful immersed coupling, and then perturb it to produce a 
nearly-successful coupling which can then be iterated to generate a successful coupling following the methods 
used to prove Theorem 13 and associated Lemmas. The immersed coupling described in Emery (2009, Lemma 
5) bears a strong family resemblance to the reflection / synchronized coupling given above (Example 3): in 
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preparation for construction of the equi-filtration coupling we first sketch a description of Emery's immersed 
coupling for BKR diffusions using the terminology of our paper. 

Application of the Tanaka formula for Brownian local time to (18) shows that |X| + |Y| = L X;(0) + L Y;(0) 
increases as the sum of the local times accumulated by X and Y at 0. In the case when (X, Y) does not begin 
at the origin, so that h = |(|Xo| + |Yo|) > 0, we can determine a real Brownian motion which drives the BKR 
diffusion by first defining a binary cadlag switch process K, which takes values or 1, changing only according 
to the following rules: 

1. K takes value on entry to the region |X| < h; 

2. K takes value 1 on entry to the region |Y| < h; 

and otherwise K is time-constant. The initial value of K is given by 

jo if|Xo|<h, 
1 otherwise. 
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Figure 4: Construction of a BKR diffusion (X, Y) from a real Brownian motion is based on a binary cadlag process 
K, which switches when it enters specific regions with boundaries determined by x,y = ±h for h. = |([X | + |Y |). 
The initial value (X ,Y ) lies on the boundary of the diamond-shaped region, and by construction the diffusion 
will never enter the interior of this region. 



We can then define a real- valued Brownian motion A by 

dA = Ksgn(Y)dX- (1 -K)sgn(X)dY. (19) 

Note that the definition of K implies that Y never vanishes when K = 1, and X never vanishes when K = 0. 
This allows us to use A to construct (X, Y) as follows: 

dX = Ksgn(Y)dA- (1 - K) sgn(X) sgn(Y) d Y, 

dY = -Ksgn(X)sgn(Y)dX-(l-K)sgn(X)dA. (20) 

Thus if K = 1 then we can construct X in terms of A (since Y does not then change sign) and then Y in terms 
of X, and conversely if K = then we can construct Y in terms of A and then X in terms of Y. In particular it 
is convenient to note the following relationships between X, Y and the driving Brownian motion A: 

dX = sgn(Y)dA, (21) 
dY = -sgn(X)dA. (22) 
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We can now sketch the construction of an immersed coupling with another BKR diffusion which also does 
not begin at the origin. Let (X,Y) be the coupled BKR diffusion, based on h = |(|Xol + |Yo|) > 0, with K 
and A defined in direct analogy to K and A. We define the coupling (a variant of the reflection /synchronized 
coupling described in Definition 3) by 

Definition 15 (Variant reflection / synchronized coupling). With the above notation, two BKR diffusions 
(X, Y) and (X, Y) (adapted to the same filtration) are said to be in reflection / synchronized coupling if their 
driving Brownian motions A and A are related by 

dA = sgn(Y) J sgn(Y) d A, (23) 

where J = — 1 until X and X first meet, after which we set J = +1. 

Note that this variant coupling does not treat X and Y symmetrically. As discussed in detail by Emery 
(2009), the variant reflection / synchronized coupling can be immersed and is almost surely successful. Given 
(X, Y), in order to reconstruct (X, Y) it is necessary to augment the filtration using an appropriate independent 
sequence of equiprobable ±1 random variables; this is because it is not possible to obtain strong solutions to 
the stochastic differential system given by (23), (20), and the analogue of (20) giving (X,Y) in terms of A. 
Note however that we can construct the paths of X and |Y| from the path of (X, Y): indeed 

Lemma 16. If (X, Y) and [X, Y) are connected by the variant reflection / synchronized coupling specified 
in Definition 15 then X and |Y| are both adapted to the natural filtration of X, and Ti = inf{t : X t = X t }, 
and the coupling time T couple — T3 = inf{t > Ti : |Y t | = |Y t | = 0} are stopping times for this filtration. In 
particular T C oupie is almost surely finite. 

Proof. For the case of X, observe that Ti is the first time t that Xo + Xo — X t hits §(Xo + Xo), and therefore 
Ti is a stopping time for the natural filtration of X. Moreover 

X = (Xo + X — XtATi ) + (XtVTi — Xtj ) ; 

hence X is also adapted to the natural filtration of X. 
For the case of |Y|, observe that |Y| satisfies 

d|Y| = sgn(Y)dY + dl_ Y;(0) , 

where L Y; ' ' is the local time accumulated by Y at 0, and we can use the Levy transform to show that it suffices 
to establish that sgn(Y) d Y is measurable with respect to the natural filtration of X. But we can employ the 
stochastic differential equation (23) determining the coupling, together with (22); 

sgn(Y)dY = -sgn(X)sgn(Y)dA 

= — sgn(X) sgn(Y) sgn(Y) J sgn(Y) d A = — sgn(X) J sgn(Y) d A . 

Thus we can deduce that J* sgn(Y) dY is adapted to the natural filtration of X, since dX = sgn(Y) d A by (21), 
and we have already shown that X is so adapted. 
Moreover after Ti we have J = 1 and X = X, hence 

d|Y|-d|Y| = sgn(Y)dY-sgn(Y)dY + dL Y:,0) -dL Y:(0) 

= - sgn(X) sgn(Y) d A - sgn(Y) d Y + d L Y:(0) - d L Y;(0) 

= -sgn(X)sgn(Y)dA-sgn(Y)dY + dL Y:(0) -dL Y:(0) = d L Y:(0) - d L Y;(0) , 

where we use (22) to cancel the two Brownian stochastic differentials. Thus after Ti it is the case that either 
|Y| |Y| for all time, or |Y| ^ |Y| for all time, and coupling occurs at the first time T 3 after Ti that |Y| and |Y| 
simultaneously hit 0. 

Now |Y| hits zero when |X| first hits |X | + |Y |, or at subsequent times when |X| attains its running supremum, 
while we have already shown that |Y| is adapted to the natural filtration of X, therefore T3 is a stopping time 
for this filtration. 

Finally, almost sure finiteness of T coup i e follows, since it is the (dependent) sum of two almost surely finite 
Brownian hitting times Ti and T3 — Ti . □ 



21 



We shall use this partial reconstruction when analyzing the equi-filtration coupling described below. 



4.2 An equi-filtration coupling for BKR diffusions 

In order to construct an equi-filtration coupling for BKR diffusions (X, Y) and (X, Y), such that neither BKR 
diffusion starts at the origin, we adopt the strategy of Section 3. Given a delayed time-change <r(t) defined in 
terms of a positive non-increasing function i|) as in (14), we can define a new driving Brownian motion A in 
terms of A via 

dA t = sgn(9 a{t) ) J t sgn(Y ff(t) )dA. (24) 

Here X (and Y) are defined in terms of A using (20); we note that J is the immersed control given in the previous 
subsection, constructed in terms of (X, Y) by setting J = — 1 till the time Ti when X first hits |(X + X ), and 
then setting J = +1; finally, Y (and X) are defined in terms of A using the analogue of (20). The use of the 
delay cr means that the system of these stochastic differential equations has a unique strong solution so long 
as neither BKR diffusion is begun at the origin. Note that there are issues in finding strong solutions to (20) 
together with the switching processes K and K if either or both of the BKR diffusions start at the origin. 

Lemma 17. Suppose (X, Y) is a BKR diffusion defined in terms of a BKR diffusion (X, Y) using the 
time-delayed stochastic differential equation (24) and the analogues of the defining equation (20) together 
with switching processes K and K, so that 

dX = sgn(Y)dA, (25) 
dY = -sgn(X)dA. (26) 

// neither BKR diffusion is begun at the origin then the resulting coupling is equi-filtration. 

Note that this definition is not symmetrical in (X, Y) and (X,Y), since the coupling control J is defined in 
terms of (X, Y). Note also that we do not assert that the coupling is successful! 

Proof. It follows from construction that (X, Y) is immersed in the filtration of (X, Y). On the other hand we 
can argue as in Lemma 11 that the reverse also holds, and hence that this coupling is equi-filtration. As in 
Lemma 11 of the argument for the case of Brownian motion with local time, the key point is to argue first 
that the trajectory of (X,Y) up to the time Ti (while J = —1) is immersed in the filtration of (X,Y), and then 
to argue that the subsequent construction (while J = +1) is also immersed. The crucial point is that Ti is a 
stopping time for the filtration of (X,Y), as noted in Lemma 16. □ 

Because of Lemma 16, it makes sense to discuss Xy 3 and |Yj 3 | defined in terms of X and Y, and in particular 
to consider the extent to which Xy 3 and |Yt 3 | differ from Xy 3 and |Yy 3 |. Moreover |Y|y 3 = Yy = 0, so control of 



|Y T ,|-|Y 



T 3 



corresponds directly to control of | Yj 3 — Yj 3 1 itself. 



Lemma 18. Suppose (X, Y) is a BKR diffusion defined in terms of a BKR diffusion (X,Y) using the time- 
delayed stochastic differential equation (24) and the analogues of the defining equations (21) and (22). 
For any 5 > 0, we can choose e € (0, h) sufficiently small so that if the time-delay cr(t) = t — (i|>(t) A t) 
is defined via i|>(t) = £ 3 /((t — e + l) 3 ) for t ^ e then 

P[|X T3 -X T3 | + |Y T3 -Y T3 |>6] < 5 (27) 

Proof. Consider the stochastic differential equation for X: 

dX = sgn(Y)dA = jBgn(Ya)Bgn(Y ff )sgn(YldA. 

Since dX = sgn(Y) d A, and since Y m and Y ss Yo-, it follows that the coupling between X and X approxim- 
ates a reflection coupling up to time Ti, and after that approximates a synchronized coupling. Calculating as 
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in Lemma 12, but recalling from Lemma 16 that X = (Xo + Xo — XtATi) + (X-tvTi — XyJ is actually adapted 
to the filtration of (X, Y), 



E 



sup <MX, X, 



s$ 4 

sC 32 I 
Jo 



E 



(J sgn(Y ff ) sgnfo) sgn(Y) - J sgn(Y))' 



dt 



[sgn(Y ff ) ^sgn(Y)]dt + 32 



[sgn(Y ff ) ^sgn(Y)]dt 



105.557... x e. 



Thus we can control the extent to which the approximate reflection coupling of X and X will deviate from the 
reflection coupling of X and X: using the Markov inequality, for any 5 > if £ < 5 3 /(12 x 105.557 . . .) then 



sup | Xt-Xt | > 5/2 



6/3. 



(28) 



We now need to control the approximation of |Y| by |Y|. We first use the almost-sure finiteness of the 
stopping time T 3 to select a constant time t max such that 

P [T 3 > t max ] < 6/3 . (29) 

It suffices to control the approximation of JgSgn(Y) dY by JgSgn(Y) dY when ^ t ^ t max . Observe that 

sgn(Y)dY-sgn(Y)dY = -sgn(Y) sgn(X) d A + sgn(Y) sgn(X) d A 

= -J sgn(Y) sgn(X) sgn(Y a ) sgn(Y a ) d A + J sgn(Y) sgn(X) d A . 

Hence (using the same definition of \\>) 
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sup < 
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sgn(Y)dY 



sgn(Y)dY) 
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s$ 48 
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[sgn(Y) + 1 sgn(Y a ] dt + 48 



[sgn(Y) #sgn(Y CT ]dt- 



+ 48 



sgn(X) ± sgn(X) 



dt 



158.336... x £ + 48 



sgn(X) ^ sgn(X) 



dt. 
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It follows that if we choose first r| then £ small enough that 
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11 < 



£ < 



£ < 



TX 



12 x 48 V 8t 



TX 
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12 x 48 x 105.557 . . . x t max 
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then 



IE 


sup < 


[fl 









sgn(Y)dY- 

o 



sgn(Y)dY) 

o 



105.557. . . x t max 2 2 /8t n 

£ + -V + \ 



158.336... x £ + 48 | - 2 ' - £ + ^T] 2 + \/ — g^T] ) < 5/3. (30) 



The lemma now follows from the inequalities (28), (29), and (30). □ 

We can now state and prove the main result of this section, that BKR diffusions begun at non-zero points 
can be coupled in an equi- filtration manner. 

Theorem 19. Given two BKR diffusions begun at different non-zero initial points, they can be coupled 
in a mutually immersed manner which succeeds in almost surely finite time, using an infinite sequence 
of increasingly rapid equi-filtration couplings, each of which approximates the variant reflection / syn- 
chronized coupling but with delays built into the reflection and synchronization couplings to render them 
equi-filtration. 

Proof. From Lemma 18 it follows that at time T 3 we can bring (X, Y) and (X, Y) arbitrarily close together with 
probability arbitrarily close to 1. Moreover Yj 3 = Yy 3 = 0, so that Yj 3 will be arbitrarily close to zero with 
probability arbitrarily close to 1. 

Restarting at time T 3 , the evolutions of (X, |Y|) and (X, |Y)| can be related to the behaviour of Brownian 
motion and its local time at 0. Specifically, (|Y|, X — jYj) (respectively (|Y[, X — |Y|)) has the stochastic dynamics 
of the absolute value of Brownian motion together with its local time at 0, at least until X (respectively X) 
hits zero. 

We can therefore apply the iterative coupling technique of Section 3 to achieve exact coupling of (X, Y) and 
(X, Y); the localization supplied by Corollary 14 implies that there is a positive probability of achieving this 
coupling before either X or X hit 0, with a uniform positive lower bound on the probability which tends to 1 
as the restarted values at T3 of X — X, Y and Y tend to zero. In the event of default (i.e. the initial variant 
reflection / synchronized coupling fails to achieve approximate coupling at T3, or X or X hits subsequent to 
T3), then the whole coupling procedure can be restarted; the initial delayed variant reflection / synchronized 
coupling can be arranged to deliver approximate coupling to an arbitrarily small tolerance with probability 
arbitrarily close to 1, and then the subsequent iterative coupling, will also have success probability arbitrarily 
close to 1. Thus it is possible to arrange that the coupling procedure will only need to be restarted a finite 
number of times before coupling is achieved. □ 



5 Conclusion 

This paper makes a detailed study of the reflection / synchronized coupling for Brownian motion together with 
local time. Using the Levy transform, relating the absolute value of Brownian motion with local time to a new 
Brownian motion together with its running supremum, it is shown that this coupling is optimal and indeed 
strictly optimal amongst all immersed couplings of Brownian motion with local time; on the other hand an 
exact calculation of the moment-generating function of the coupling time (for the non-singular case when both 
reflection and synchronized stages are required) indicates that the coupling is not maximal. The reflection 
/ synchronized coupling cannot be an equi-filtration coupling, for reasons associated with Tanaka's famous 
example of a stochastic differential equation which has weak solutions but no strong solutions. However a 
simple modification, involving deterministic time-delays, delivers a coupling composed of an infinite sequence 
of increasingly rapid time-delayed approximate couplings. This coupling is an equi-filtration coupling, and 
almost surely succeeds in finite time. 

It would be an interesting and significant exercise to describe an exact construction of the maximal coupling, 
based on a continuous-time modification of Pitman (1976)'s construction for discrete-time discrete-state-space 
Markov chains. One might expect, as a byproduct, a rigorous (rather than merely numerical) proof of non- 
maximality for the reflection / synchronized coupling. We leave this for a further investigation, as well as the 
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question of whether in the singular case (a) the truncated reflection / synchronization coupling continues to 
be optimal amongst immersed couplings, and (b) continues not to be maximal (we believe both of these to 
be the case). Rather than pursuing these questions, the paper describes an application of this detailed study 
of coupling of Brownian motion with local time. Emery (2009) raised the question of whether it would be 
possible to make an explicit construction of an equi-filtration coupling of the BKR diffusion (Benes et al., 
1991), specifically when neither of the two copies to be coupled are started at the origin. The deterministic 
time-delay trick is applied to Emery (2009) 's immersed coupling of the BKR diffusion, the variant reflection / 
synchronized coupling, to produce an almost surely successful coupling which works by iteration of approximate 
equi-filtration couplings. The strategy closely follows the case of Brownian motion with local time, but requires 
somewhat more detailed analysis. 

Further questions include the following: 

1. Is Emery (2009) 's immersed BKR coupling optimal amongst all immersed couplings of BKR diffusions? 
We would be surprised if this were the case, since there are two distinct variants of Emery (2009) 's 
coupling depending on whether the construction is based on the X component or the Y component: it 
seems implausible that optimal couplings would permit such a symmetry. It is of course possible that 
there is no optimal immersed coupling: it may be the case that some strategies work better for short 
time while others work better for long time. 

2. Reverting to the case of Brownian motion with local time, is it possible to couple Brownian motion 
together with local times accumulated at two or more distinct points? This seems to be a hard question. 
Analogous generalizations have been carried out for the case of coupling Brownian motions together with 
iterated time-integrals (Kendall and Price, 2004); however there is a useful linearity in the time- integral 
case which is not present here. The question of coupling finite sets of local times could be viewed as a 
question of whether one could couple a finite version of the celebrated Brownian burglar (Warren and 
Yor, 1998; Aldous, 1998). 

3. There is interesting territory to be explored in the realm of couplings which fall short of being maximal 
but yet are not immersed. One example in applied stochastic process theory is supplied by Smith (2011), 
who investigates the mixing time of a simple Gibbs sampler on the unit simplex using a two-stage 
coupling of which the first is immersed (Markovian, in Smith's chosen terminology) while the second 
couples an associated partition process anticipatively. This non-immersed coupling allows Smith to give 
an affirmative answer to a conjecture by Aldous concerning the mixing time of this Gibbs sampler. 
Arguably Sigman (2011)'s perfect simulation algorithm for super-stable M/G/c queues can be put in 
the same category, as this depends on coupling service times not according to time of arrival of customer 
but according to time of start of service. It would be likely to be most illuminating if one could discover 
simple Brownian coupling problems for which gains of a similar kind can be made. 

4. It seems clear that the techniques of this paper can be generalized to show that immersion couplings 
of suitably regular diffusions can always be approximated by equi-filtration couplings, and it would be 
interesting to see a fully rigorous proof in a case where the qualification "suitably regular" is given a 
pleasant and natural meaning. 

The underlying reflection / synchronization coupling for Brownian motion together with local time is 
extremely simple, and lends itself to rather complete calculation. Not only is it an example of the general 
programme of coupling Brownian motion together with appropriate functionals (Ben Arous et al., 1995; Kendall 
and Price, 2004; Kendall, 2007, 2010), but also it can be viewed as a basic coupling strategy that, like the 
reflection coupling of Brownian motion (Lindvall, 1982) has the potential to serve as a model in much more 
general situations. The application to the BKR diffusion in this paper illustrates this point; it is hoped that 
the calculations described above will facilitate the use of the reflection / synchronization coupling as a building 
block in other applications of coupling to probability theory. 



25 



Acknowledgements 



I gratefully acknowledge the contribution of Prof. Michel Emery, whose question concerning BKR diffusions 
encouraged me to develop the initial idea of coupling Brownian motion together with local time, and who 
commented insightfully on an early draft of this paper. 

References 

Aldous, D. J. (1998). Brownian excursion conditioned on its local time. Elect. Comm. in Probab 3, 79-90. 

Beghdadi-Sakrani, S. and M. Emery (1999). On certain probabilities equivalent to coin-tossing, d'apres Schach- 
ermayer. Seminaire de Probabilites XXXIII 33, 240-256. 

Ben Arous, G., M. Cranston, and W. S. Kendall (1995). Coupling constructions for hypoelliptic diffusions: 
Two examples. In M. Cranston and M. Pinsky (Eds.), Proceedings of Symposia in Pure Mathematics, 
Volume 57, Providence, RI Providence, pp. 193-212. American Mathematical Society. 

Benes, V. E., I. Karatzas, and R. W. Rishel (1991). The separation principle for a Bayesian adaptive control 
problem with no strict-sense optimal law. In Applied Stochastic Analysis, Stochastics Monographs, pp. 
121-156. New York: Gordon & Breach. 

Burdzy, K. and W. S. Kendall (2000, May). Efficient Markovian couplings: examples and counterexamples. 
The Annals of Applied Probability 10(2), 362-409. 

Chen, M.-F. and S. F. Li (1989). Coupling methods for multidimensional diffusion processes. The Annals of 
Probability 17(1), 151-177. 

Connor, S. (2009, August). Optimal co-adapted coupling for a random walk on the hyper-complete-graph. 
arXiv 0908.2038, 16 pp. 

Connor, S. B. (2007). Coupling: Cutoffs, CFTP and Tameness. Phd thesis, University of Warwick. 

Connor, S. B. and S. D. Jacka (2008). Optimal co-adapted coupling for the symmetric random walk on the 
hypercube. Journal of Applied Probability ^5(1), 703-713. 

Doeblin, W. (1938). Expose de la Theorie des Chaines simples constants de Markoff a un nombre fini d'Etats. 
Revue Math, de I'Union Interbalkanique 2, 77-105. 

Emery, M. (2005). On certain almost Brownian nitrations. Ann. Inst. H. Poincare Probab. Statist. ^1(3), 
285-305. 

Emery, M. (2009). Recognizing Whether a Filtration is Brownian: a Case Study. In C. Donati-Martin, 
M. Emery, A. Rouault, and C. Strieker (Eds.), Seminaire de Probabilites XLII, Volume 1979 of Lecture 
Notes in Mathematics, pp. 383-396. Berlin, Heidelberg: Springer Berlin Heidelberg. 

Goldstein, S. (1978). Maximal coupling. Zeitschrift fur Wahrscheinlichkeitstheorie und Verwe Gebi- 
ete 46(2), 193-204. 

Griffeath, D. (1975). A maximal coupling for Markov chains. Zeitschrift fur Wahrscheinlichkeitstheorie und 
Verwe Gebiete 31, 95-106. 

Ikeda, N. and S. Watanabe (1981). Stochastic differential equations and diffusion processes. Amsterdam: 
North Holland / Kodansha. 

Ito, K. (1975). Stochastic differentials. Applied Mathematics and Optimization 1, 374-381. 

Kendall, W. S. (1998). From Stochastic Parallel Transport to Harmonic Maps. In Juergen Jost and Wilfrid 
S Kendall and Umberto Mosco and Michael Roeckner and Karl-Theodor Sturm (Ed.), New Directions in 
Dirichlet Forms, pp. 49-115. Providence, RI Providence: American Mathematical Society. 



26 



Kendall, W. S. (2001). Symbolic Ito calculus: an ongoing story. Statistics and Computing 11, 25-35. 

Kendall, W. S. (2007, May). Coupling all the Levy stochastic areas of multidimensional Brownian motion. 
The Annals of Probability 35(3), 935-953. 

Kendall, W. S. (2009, September). Brownian couplings, convexity, and shy-ness. Electronic Communications 
in Probability 14 (Paper 7), 66-80. 

Kendall, W. S. (2010). Coupling time distribution asymptotics for some couplings of the Levy stochastic area. 
In N. H. Bingham and C. M. Goldie (Eds.), Probability and Mathematical Genetics: Papers in Honour 
of Sir John Kingman, Chapter 19, pp. 446-463. Cambridge: Cambridge University Press. 

Kendall, W. S. and C. J. Price (2004). Coupling iterated Kolmogorov diffusions. Electronic Journal of 
Probability 9(Paper 13), 382-410. 

Kuwada, K. (2009). Characterization of maximal Markovian couplings for diffusion processes. Electronic 
Journal of Probability 14, 633-662. 

Kuwada, K. and K.-T. Sturm (2007). A counterexample for the optimality of Kendall-Cranston coupling. 
Electronic Communications in Probability 12, 66-72. 

Lindvall, T. (1982). On Coupling of Brownian Motions. Technical report 1982:23, Department of Mathematics, 
Chalmers University of Technology and University of Goteborg. 

Lindvall, T. (1991). W. Doeblin 1915-1940. The Annals of Probability 19(3), 929-934. 

Lindvall, T. (1992). Lectures on the coupling method. Wiley Series in Probability and Mathematical Statistics: 
Probability and Mathematical Statistics. New York: John Wiley & Sons Inc. 

Lindvall, T. and L. C. G. Rogers (1986). Coupling of multidimensional diffusions by reflection. The Annals 
of Probability 14(3), 860-872. 

Pitman, J. W. (1976). On coupling of Markov chains. Zeitschrift fur Wahrscheinlichkeitstheorie und Verwe 
Gebiete 55(4), 315-322. 

Revuz, D. and M. Yor (1991). Continuous martingales and Brownian motion, Volume 293 of Grundlehren 
der Mathematischen Wissenschaften [Fundamental Principles of Mathematical Sciences]. Berlin: 
Springer-Verlag. 

Rogers, L. C. G. (1999). Fastest coupling of random walks. The Journal of the London Mathematical 
Society (Second Series) 60(2), 630-640. 

Rogers, L. C. G. and D. Williams (1987). Diffusions, Markov processes, and martingales, Volume II. 
Chichester / New York: John Wiley & Sons Ltd. 

Rosenthal, J. S. (1997). Faithful Couplings of Markov Chains: Now Equals Forever. Advances in Applied 
Mathematics 18(3), 372-381. 

Sigman, K. (2011). Exact Simulation of the Stationary Distribution of the FIFO M/G/c Queue. Journal of 
Applied Probability 48A, 209-213. 

Smith, A. (2011, July). A Gibbs Sampler on the n-Simplex. arxiv preprint arxiv: 1 107, 16. 

Sverchkov, M. Y. and S. N. Smirnov (1990). Maximal coupling for processes in D[0,oo]. Dokl. Akad. Nauk 
SSSR 311(b), 1059-1061. 

Thorisson, H. (1994). Shift-coupling in continuous time. Probability Theory and Related Fields 99(A), 
477-483. 

Thorisson, H. (2000). Coupling, stationarity , and regeneration. New York: Springer-Verlag. 



27 



Warren, J. and M. Yor (1998). The Brownian burglar: conditioning Brownian motion by its local time process. 
Semmaire de Probabihtes XXXII 32, 328-342. 



Department of Statistics, University of Warwick, Coventry CV4 7AL, UK 
Email address: w.s.kendall@warwick.ac.uk 



28 



